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Leggett-Garg inequalities (LGI) test the correlations of a single system measured at different times. 
Violation of LGI implies either the absence of a realistic description of the system or the impossibility of 
measuring the system without disturbing it. We investigate the violation of the Leggett-Garg inequality for a 
two level system under decoherence in a non-Markovian dephasing environment. We discuss the 
non-Markovian dynamics of the violation of LGI at zero temperature and also at finite temperature for 
different structured environments. An enhanced quantum coherence is shown through the violation of 
Leggett-Garg inequahty in the strong non-Markovian regime of the environment. 

There has been recent attention and renewed interests in the investigation of Leggett-Garg (LG) inequalities. 
The original motivation' for these inequalities was to test quantum coherence in macroscopic systems. LG 
inequalities also play a useful role for microscopic systems as an indicator of nonclassicality due to the 
existence of superposition states and their collapse under measurement in the quantum regime. The LG inequal- 
ities were derived under following two assumptions (Al) macrorealism: performance of a measurement on a 
macroscopic system reveals a well-defined pre-existing value {A2) noninvasive measurability: this pre-existing 
value can be measured without disturbing the system. Violation of Leggett-Garg inequalities (LGI) implies either 
the absence of a realistic description of the system or the impossibility of measuring the system without disturbing 
it. Quantum mechanics (QM) violates the inequalities on both accounts. In recent years, theoretical study of QM 
violation of LGI has been carried out for a variety of systems like electron transport through quantum dot (to 
distinguish classical and quantum transport)^, for investigating nonclassicality in opto-electro-mechanical sys- 
tems', to investigate quantum coherence in oscillating neutral kaons* and in biological light-harvesting protein 
complexl The first measured violation of LGI was announced'' by Palacios-Laloy et al. and within a few years, 
violations had been probed in a diverse range of physical systems such as photons' defect centres in diamond'", 
nuclear magnetic resonance" '-, phosphorus impurities in silicon'^. 

Violation of LGI is associated to the presence of coherence in the quantum system which inevitably 
interacts with its surrounding environment (reservoirs or baths) and the coherence is lost in general due 
to this unavoidable interaction. Generally, the non-unitary evolution of the reduced-density matrix of the 
system is obtained after taking partial trace of the total system-plus-reservoir density operator pj{t) over the 
reservoir degrees of freedom. Some approximations are often made in the derivation of master equation for 
the system's reduced density matrix. The most important approximations'''"' are the weak coupling or Born 
approximations, assuming that the coupling between the system and the reservoir is small enough to justify a 
perturbative approach, and the Markov approximation, assuming that the correlation time of the reservoir is 
very short compared to the typical system response time so that the reservoir correlation function is assumed 
to be delta correlated in time. Although, the use of the Markovian approximation is justified in a large variety 
of quantum optical experiments, one should notice that non-Markovian effects are crucial, for example, for 
high-speed quantum communication where the characteristic time of the relevant system becomes comparable 
with the reservoir correlation time, or if the environment is structured with a particular spectral density, for 
example, for quantum systems embedded in solid-state devices, where memory effects are typically non- 
negligible. In these cases, the dynamics can be substantially different from the Markovian dynamics. Due 
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to their fundamental importance in quantum information proces- 
sing and quantum computation, non-Markovian quantum deco- 
herence have attracted much attention in recent years"* one of 
the main purposes of which in the long run is to engineer different 
types of (artificial) reservoirs and couple them to the system in a 
controlled way. 

The evolution equation of an open system density matrix is 
governed by the reduced Liouville equation or called the quantum 
master equation that can be Markovian or non-Markovian. In the 
Markovian case, an useful procedure to calculate the two-time 
(multi-time) correlation functions (CPs) for open (dissipative) 
quantum systems is the so-called quantum regression theorem 
(QRT)" '' that gives a direct relation between the time evolution 
equation of the single-time expectation values and that of their 
corresponding two-time (multi-time) CPs. So, knowing the master 
equation and the time evolution of the reduced density matrix, one 
is able to calculate the single-time expectation values as well as the 
two-time (multi-time) CPs using QRT. But knowing the time 
evolution of the reduced density matrix is not sufficient to calculate 
the temporal (two-time or multi-time) correlation function (CPs) 
of system observables in non-Markovian case, for non-Markovian 
open (dissipative) quantum systems, the QRT is not valid in gen- 
eraP'^*". Leggett-Garg inequalities test the correlations of a single 
system measured at different times, so one needs to calculate these 
two-time correlation functions "{Q{tj)Q(tj))" of an observable Q 
for this open quantum system. Two-time (multi-time) correlation 
functions (CPs) of an open quantum system in itself are important 
physical quantities A very useful evolution equation that allows 
systematically calculating the two-time correlation functions (CPs) 
of system operators for non-Markovian open (dissipative) 
quantum systems is derived*" recently. The derivation is based 
on perturbative quantum master equation approach, valid to sec- 
ond order in the system- environment interaction. We wUl use this 
method to calculate the two-time correlation functions for a two 
level system in a structured non-Markovian dephasing envir- 
onment, and finally investigate the non-Markovian dynamics of 
LGI-violation (coherence) for this specific open quantum system. 
The two-time CPs obtained using the derived evolution equa- 
tion'"'''" in the weak system-environment coupling case for this 
non-Markovian pure-dephasing model happen to be the same as 
those obtained from the exact evaluation. However, the two-time 
correlation functions so obtained significantly differ from the non- 
Markovian two-time CPs obtained by wrongly directly applying 
the quantum regression theorem (QRT), a useful procedure to 
calculate the two-time CPs for weak-coupling Markovian open 
systems. 

The simplest LGI can be constructed as follows. Let Q{t) be the 
observable of a two level system such that, whenever measured, it is 
found to take a value + 1 or — 1, depending on whether the system is 
in state | + ) or | — ). We then perform three set of experimental runs 
such that in the first set of runs Q is measured at times ti and t2 = f i + 
Af; in the second, at ti and ti = ti + 2At; in the third at t2 and (here 
{3 > t2 > fi). Prom such measurements, it is intuitive to determine 
the two-time correlation functions Kjj = (Q(f/)Q(fi)) where tj > ti. 
Leggett and Garg' adopted the standard argument leading to a Bell- 
type inequality, with times tj and tj playing the role of apparatus 
settings. Por any set of runs corresponding to the same initial state, 
any individual Q(f) has the same definite value {assumption Al), 
irrespective of the pair Q(t,)Q(f,) in which it is measured; i.e., the 
value of Q{tj) or Q(t,) in any pair does not depend {assumption A2) 
on whether any prior or subsequent measurement has been made on 
the system. Consequently the combination Q{t2)Q{ti) + Q(t3)Q(t2) 
^ Qit3)Q{ti) has an upper bound of + 1 andlowerboundof — 3. If all 
the individual product terms in this expression are replaced by their 
averages over the entire ensemble of such sets of runs, the following 
form of LGI is then obtained: 



K3=K2l-hK32-K3l<l (1) 

The same LGI can also be derived through probability argument*^. 
With similar arguments one can also derive a Leggett-Garg inequal- 
ity for measurements at four different times, ti, t2, t^ and 14 = ti + 
3 At given by 

K4 = K2i -|--fC32 -I--K43 — K41 <2 (2) 

The two-time operators Q(t/)Q(f,) are not Hermitian operators, 
and the two-time correlation functions Kjj = {Q{tj)Q{tj)) are in 
general complex quantities. We consider the symmetrised com- 
bination of the two-time correlators: Kjj = {{Q{tj), Q((,)})/2, to 
indentify them with physical expectation values of the two-time 
measurements. Please note that the symmetrised operator {{Q{tj), 
Q{tj)})/2 = {Q{tj)Q{tj) + Q{tj)Q{tj))/2 is Hermitian whose expecta- 
tion value can be associated to the average value of the two-time 
measurements (see Supplementary Information and''''). Hence, we 
will only consider the real part of K3 and K4 to investigate the 
Leggett-Garg inequalities for a two level system interacting with a 
non-Markovian dephasing environment. Recently, a general frame- 
work for understanding the influence of non-unitary evolution on 
maximal violations of the LGI was given*'' for a class of Markovian 
decoherence channels. Using LG inequalities, another proposal''^ 
was given to test quantum coherent dynamics of a two-qubit sys- 
tem undergoing through a superradiant decay under a common 
reservoir, Born-Markov master equation was used for the open 
system dynamics. Por the experimental realization, they provided 
a feasible scheme, which consists of two quantum dots coupled to 
nanowire surface plasmons. Dynamics of the violation of the 
Leggett-Garg inequality for a two level system under decoherence 
in a non-Markovian dephasing environment, to our knowledge, 
have not been presented in the literature. We wiU briefly discuss 
the pure-dephasing spin-boson model and the perturbative method 
to calculate the non-Markovian evolution equation for reduced 
density matrix and two-time CPs valid up to second order in 
system-environment coupling strength. Then we present our 
numerical results to investigate Leggett-Garg inequality under the 
non-Markovian decoherence dynamics, where we show an 
enhanced quantum coherence through the violation of Leggett- 
Garg inequality in a suitably tuned non-Markovian regime of the 
environment. 

Results 

Pirst, we discuss on the evolution equation of non-Markovian two- 
time correlation functions. We consider a general class of open 
quantum systems modeled by the total HamUtonian of the system 
plus reservoir as 

H = Hs -\- + Hi 

k k 

where Hg and Hf, are system and environment HamUtonian, respect- 
ively, and Hj describes the interaction between the system and the 
environment. Here the operator L acts on the Hilbert space of the 
system, a'j^ and O/^ are the creation and the annihilation operators on 
the bosonic environment Hilbert space, and gk and ojj^ are the coup- 
ling strength and the frequency of the fcth environmental oscillator, 
respectively. The reduced density operator of the system can be 
obtained from the density operator of the total system by tracing 
over the environmental degrees of freedom ps{t) = Tr^ {prit))- 
The total density operator is governed by the quantum evolution: 

Pj.{t) = e-r'"'pT{0)er'"'.One can obtain a time-convolutionless non- 
Markovian master equation at finite temperature, valid to second 
order in the system-environment interaction strength 
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dpM 1 r 

-IfAHs, Psit)] 



dt ih ' 



no =2 



- ' dT{a{t-T)[L^L{T-t)p,{t)-L{-z-t)p,{t)L^] 

JO 

+ jS(t-T) [LV[x-t)ps[t)-i\z-t)ps{t)L]+l\.c.}, 



where 



(5) 



(6) 



are the environment CPs with a(f — T) = y^^.g* 0^,(1) 



and ^(f-T) = ^^^^*a^(0^^,^rflt'(T)^^ , where 5^(0 = 

fliie^™'' and a|(t) =fl|^e""''' are the environment operators in the 
interaction picture, and the symbol {■ ■ ■)k denotes taking a trace with 
respect to the density matrix of the thermal reservoir (environment). 
Here ni^ = {exp{ha)f;/kBT) — l)^^ is the thermal mean occupation 
number of the bosonic environment oscillators. 

Using open quantum system technique and quantum master 
equation approach, one can also calculate the non-Markovian evolu- 
tion equation of the two-time CPs of system observables Q{tj) and 
Q{tj) (for finite temperature environment), valid to second order in 
system-environment coupling as'*° 

d{Q{tj)Q{t,))/dtj 

= (i/fi)Tr,({[H,,Q]}(f,)Q(f,)p(0)) 



dzTis 

{a* [tj - t) {r (t - tj) [Q,L] } [tj) Q(t,)p(0) 
+ ^{tj - t) { [L\Q]L{T-ti) } (t,) Q(f,)p(0) 
+ P*{tj-T){i{T-t,) [Q,L^]}{t,)Q{t,)p{0) 
+ p{tj-T){[L,Qp{,-tj)}{t,)Q{u)p{0) 



(7) 



dzTrs 

{a{t,-T){[L\Q] } {tj){[Q, i(T-f.)] }{t,)p{0) 
+ p{tj-T){[L,Q]}{t,){[Q:LHT-t,)] }(f,)p(0)). 

where L{t) = exp{iHst/h)L exp( — iHst/h) is the system operator of 
the interaction Hamiltonian (Eq. 3) in the interaction picture with 
respect to Hg and Q = Q(0). In the derivation of Eqs.(4) and (7), we 
have considered the factorized initial system-bath state p j-(O) = ps(0) 
® Pr, where pR = exp{ — PHR)/Trlexp{ — I^Hr)]. Quantum regres- 
sion theorem for this non-Markovian evolution is not applicable due 
to the presence of the last integral term in Eq.(7). For pure dephasing 
spin-boson model L = = V, and 

Hs = {hwo/2)(T,, (8) 

then the non-Markovian master equation for the reduced system 
density matrix becomes"" 

= ^W.,Ps{t)]-^[Ps{t)-c..Psit)rr.]. (9) 

The non-Markovian effect in the master equation (9) is taken into 
account by the time-dependent coefficient r{t) given by 



dT(0(eff (t- t) + attt{t - t)) , 



where a.cff(t ^ t) is the effective bath CP given by'" 

«eff (f — t) =Ol{t — z)+P{t — z) 



(10) 



(11) 



where a{t — t) and /?(f — t) are defined in Eqs.(5) and (6). The single- 
time expectation value of an operator can be obtained as 



(Q(t))=Trs8«(Q(0Pr(0)) 

= Trs(x,i^(e'^Q(0)e-'Trp^,(0)^ 

= Trs®;i(Q(O)pj.(0)=Trs(Q(O)ps(0) 



(12) 



Using Eq.(7), one can also calculate the non-Markovian evolution 
equation of the two-time CPs of system observables as 



— (ffx {t,) CTyit.) )=-r{tj) (a, {t,) (Jyit.) ) - COo{ay {tj) Gyiti) ) 

-ri{tj,ti){ay{tj)a,{t,)), 

^ (ff;, (tj) aM ) = -r{tj){ay {tj)a,{t,) ) + C0o{a, (tj) a,{t,) ) 
-ri{tj,ti){a,{tj)ay{ti)), 

(ff.(t;)ff.(f,))= -r(t;)(ff,(t;)ff,(t,))-fflo(ff;,(t;)ffx(t,)) 

+ r,{t,A){rTy{t,)ay{t.)), 

{oy (t,-) ay {t,)) = -T{t,) {ay {t,) ay{U) ) + ojo{ay (tj) ay{u) ) 
+ T,{tj,t,){a,{t,)a,{t,)), 



dt 



dt 



where 



dzaLf.if{tj - 



(13) 



(14) 



(15) 



(16) 



(17) 



To calculate the functions Y{t) and Ti{tj, t,), we need to evaluate the 
environment correlation function: 



adf{t-z)-- 



dwJ{oj){coth{ho}/2kBT)cos[o}{t — z) 
-i sin[w{t — z)]} 



where /(cu) = Ej;jg;tj^5(cij — co^-) is the spectral density of the envir- 
onment. We may consider any spectral density to characterize the 
environment, here we consider the following spectral density 

/(co)=ycoQ'"'exp(-^) (19) 

where y is the coupling strength between the system and the envir- 
onment and A is the cutoff frequency. When s < 1, 1, and >1, the 
corresponding environments are sub-Ohmic, Ohmic, and super- 
Ohmic respectively. Let us now discuss the non-Markovian charac- 
teristics of these environments through the dynamics given by 
Eq.(9). Recently important steps towards the development of a gen- 
eral consistent theory of non-Markovian quantum dynamics have 
been made and several measures have been proposed^^ which try to 
define the border between Markovian and non-Markovian quantum 
evolution and to quantify memory effects in the open system 
dynamics. It is important to mention here that non-Markovianity 
can be used as a resource in the new quantum technologies. One can 
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■ non-Markovian: A=l/3 to. 



- non-Markovian: A=3 co^ 

- non-Markovian: A=30 co„ 



(a) Sub-Ohmic 



(b) Ohmic 



(c) Super-Ohmic 




Figure 1 | We plot the trace distance D(pi(t), p2(t)) at zero temperature for three different values of s = 0.5, s = 1, s = 3 corresponding to sub-Ohmic, 
Ohmic and super-Ohmic environments with / = 0.05. For each environment, we consider three different cutoff frequencies A = coo/3, A = ojq, 
and A = 30(Uo. 



induce non-Markovianity, in the spirit of reservoir engineering, to 
improve quantum protocols such as quantum metrology and 
quantum key distribution''' Here, we consider a particular 
approach to demonstrate the non-Markovian dynamics of the two 
level system, based on the trace distance between quantum states^'*. 
Quantum memory effects can be visualized through the dynamics of 
the trace distance D(pi, pj) = (1/2) Tr \pi — P2! between two 
quantum states pi and p2. This quantity can be interpreted as a 
measure for the distinguishability of the two states. Markovian pro- 
cesses tend to continuously reduce the distinguishability of physical 
states, which means that there is a one way flow of information from 
the open system to its environment. In order to have non-Markovian 
effects, there must be, for some interval of time, an information flow 
from the environment back to the system. In view of this interpreta- 
tion the characteristic feature of a non-Markovian quantum process 
is the increase of the distinguishability, that is a reversed flow of 
information from the environment back to the open system. The 
information flowing from the environment back to the system allows 
the earlier states of the system to have an effect on the later dynamics 
of the system, that is, it allows the emergence of memory effects. The 
non-Markovian dynamics of the two-level system is illustrated in 
Fig. 1 where we plot the trace distance D{pi{t), Pzit)) for the follow- 
ing pair of initial states 



Pi(0) = |+)( + |, 



P2{0) - ^{\+) + \ -)){{ + \ + {-\), (20) 



where where | +) and | — ) are the eigenstates of CTj.. Fig. 1 shows that 
the non-Markovian memory effect is enhanced for suitable choices of 
the cutoff frequencies for Sub-Ohmic and Ohmic environment, 
whereas the strong memory effect under a super-Ohmic envir- 
onment is less sensitive to the cutoff frequency. We plot in Fig. 2, 
the noise power spectrum S(co) for different type of environments. 
S(a)) is just the Fourier spectrum of {(Tx(t + T)o"-t(0) calculated using 
Eqs.(13-16). We compare the non-Markovian noise power spectrum 
S(m) with its Markovian counterpart calculated by using the 
quantum regression theorem (QRT). The coherent peaks of the 
Fourier spectrum are higher and widths are narrower in the non- 
Markovian evolution. It is important to mention here that the noise 
power spectrum S{co) was directly measured in the first experimental 
test of LGI violation'', and two-time correlation functions were then 
computed by inverse Fourier transformation of S{oj). Next, we 
investigate the dynamics of LGI for a two-level system under deco- 



herence in a non-Markovian dephasing environment with 0 = 0^ 
and Kjj = {cixit^aj^t^) using Eq.(7). We stress that the non- 
Markovian effects are crucial for open quantum systems interacting 
with an environment structured with a particular spectral density. 
We show that the violation of LGI, or in other words the quantum 
coherence dynamics can be controlled with suitable choices of the 
environment spectral density. The initial environment state is con- 
sidered to be in the thermal equilibrium state and the system is 

arbitrarily chosen to = ^ | -h ) -h - 1 - ), hence pM = |^) <*| ■ 

Then we calculate K3 and K4 in three different cases with different 
structured environments characterized by the spectral density (19). 
In Fig. 3, we plot the real part of and for different values of the 
environmental parameter s. We show the dynamics of this inequality 
for different structured environments by varying the cutoff fre- 
quency A. In Figs. 3(a) and (d), we consider a sub-Ohmic (s = 
0.5) zero temperature environment. In Figs. 3(b) and (e), we consider 
vacuum reservoir with Ohmic (s = 1) spectrum. Figs. 3(c) and (f) 
show the dynamics of and for a super-Ohmic (s = 3) vacuum 
bath. In general, the system dynamics goes beyond classical descrip- 
tion (violation of LGI) for short measurement intervals At. From 
Figs. 3(a) and (d), we notice that the quantities and K4 are oscil- 
latory and violate the LG inequalities with the strongest violation 

Tt 

occurring at the first maxima, at At = , with cut-off frequency A 

3oJo 

= cOo/3. The maximum QM value of K3 or K4 depend on the value of 
the cut-off frequency A of the spectral density. We observe an 
enhanced dynamics of quantum violation of LGI {K^ > 1 and 
> 2) at lower cut-off frequencies for sub-Ohmic and Ohmic envir- 
onment. The enhanced quantum coherence dynamics is due to the 
stronger non-Markovian memory effect when the cut-off frequency 
is small. For sub-Ohmic and Ohmic spectrum, we see a reduced 
violation of LGI (Figs. 3(a), (d) and (b), (e)) as we increase the 
cut-off frequency. Fig. 3(c) and (f) indicate that the LGI violation 
is less sensitive to cutoff frequency in case of super-Ohmic (5 = 3) 
environment, and we have wide range of At where LGI is violated in 
case of super-Ohmic spectrum. This is also clear from Fig. 1, where 
we see an enhanced non-Markovian memory effect for Sub-Ohmic 
and Ohmic environment when the cut-off frequency is small, also 
the trace distance D(pi(t), pzit)) is less sensitive to the cutoff fre- 
quency in case of super-Ohmic environment. Next, we go to the 
finite-temperature case for which h;i = {exp{hcOk/kBT) — l)^^ j^Q. 
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(a) Sub-Ohmic 



■ Markov: QRT 

■ non-Markovian 

(b) Ohmic 



(c) Super-Ohmic 



^3 





-10 -5 0 5 10 

UJ 

Figure 2 | We plot the noise spectrum S(oj) at zero temperature for three different values of s = 0.5, s = 1, s = 3 corresponding to sub-Ohmic, 
Ohmic and super-Ohmic environments with cutoff frequency A = g>o/3, coupling strengths y = 0.05 and f = 0.1. We compare the non-Markovian noise 
power spectrum S(f/j) with its Markovian counterpart calculated by using the quantum regression theorem (QRT). 

In Fig. 4, we show the effect of temperature on LG inequality for the three different values of the environmental parameters "s" with fixed 
system arbitrarily chosen in an initial state |T) and the environment value of A = cuq. We find that at low temperature {k^T = O.OltOo), 
initially in the thermal equilibrium state. We calculate and K4 for revival of the violation of LGI occur, but at higher temperature (fc^T 



■ non-Markovian: A=l/3 o), 



■ non-Markovian: A=3 co^ 

■ non-Markovian: A=30 ro, 




ID 





(d) Sub 


-Olimic 









1 2 
At 




Figure 3 | We plot real parts of K3 and K4 at zero temperature for three different values of s = 0.5, s = 1, s = 3 corresponding to sub-Ohmic, 
Ohmic and super-Ohmic environments with 7 = 0.05. For each environment, we consider three different cutoff frequencies A = coq/B, A = coq, and 
A = 30coo. 
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Figure 4 | We plot real part of K4 at finite temperature for three different values of s = 0.5, s = 1, s = 3 corresponding to sub-Ohmic, Ohmic and super- 
Ohmic environments with / = 0.05 and A = oJo. For each environment, we consider three different temperatures k^T = O.Olcoo, k^T = ojq, and 
kgT = Scoq. 



= 5ftJo) we see only one-peak violating the LGI (for very short mea- 
surement interval Af), and the system generally goes to classical 
regime at higher temperatures. Fig. 5 describes an explicit comparison 
between the LGI violation for the open system dynamics and the close 
dynamics (in absence of the environment) under the Hamiltonian 
(8). Two-time correlation functions for the closed system are calcu- 
lated following the calculation in the supplementary material. We see 
that the LGI violation is reduced when the coupling strength y of the 
two level system with its environment is gradually increased. 



Discussion 

Non-Markovian quantum decoherence dynamics have attracted 
much attention in recent years"* with interesting memory effects, 
e.g., non-divisibility of dynamical maps, back flow of information, 
dissipative dynamics with negative decay rate, non-monotonous 
increase of entropy of an open system, breakdown of quantum 
regression theorem etc., they all have fundamental importance in 
quantum information processing and quantum computation. We 
investigate the dynamics of quantum coherence through LG inequal- 



(a)Sub-Ohmic 




■ Closed system 
■7^0.05 
■Y=0.2 
■y=0.5 
y=l 



(b) Ohmic 




(c)Super-Ohmic 




Figure 5 | We compare the open system LGI violations with that corresponding to close system case. We plot the real part of K4 for the closed system and 
that in presence of the reservoir at zero temperature. We consider three different values of s (s = 0.5, s = 1, s = 3) corresponding to sub-Ohmic, 
Ohmic and super-Ohmic environments with cutoff frequency A = coo/3. We consider four different values of the coupling strengths y = 0.05, 0.2, 
0.5, 1.0. 
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ities for a two-level system under decoherence in a non-Markovian 
dephasing environment where we show the dynamics of the violation 
of LGI at zero or finite temperature for structured environments with 
different spectral densities (19). We have utilized an useful method'"' 
to calculate the dynamics of two-time correlation function and the 
noise spectrum for non-Markovian open quantum systems. An 
enhanced violation of LGI is observed in the strong non- 
Markovian regime when we lower the cutoff frequency of the envir- 
onment. It is important to understand that how does the memory 
(non-Markovianity) of the environment quantitatively affect the 
violation of LGI. Generally, the non-Markovian nature of the system 
appears when the system gets correlated with the environment, so it 
is not obvious how does this system-environment correlation is 
related with the measurement correlation when the system is mea- 
sured at different times. Our aim is to quantitatively address this 
question which is missing in the Uterature. 

Methods 

The two-time correlation functions {Q{tj)Q{ti)} are calculated by solving the linear 
coupled differential equations (Eqs.13-16) with the initial values are given by fol- 
lowing single- time expectation values 

(<T,(fi)a,(t,))=Trs(a^^Ps('/)) = l. (21) 

(<7,,(t,>,,(f,)) = Trs(<7>s(t,))=l, (22) 

(<T,(t,>,,(i,)) = Trs(<7,<7,,ps((,)) = i(a,(f,)>, (23) 

(<7,(f,>,(ti)> =Trs(a,a.Ps((i)) = -i((7,(f,)), (24) 
which are calculated using Eq.(12). The trace distance is defined as 

0(Pi,P2) = (l/2)Trlpi-p2|, (25) 

where the modulus of an operator A is defined by |A| — V A''' A. 
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